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Abstract 

Nuclear binding energies are investigated in two variants of the 
Skyrme model; the first replaces the usual Skyrme term with a term 
that is sixth order in derivatives, and the second includes a potential 
that is quartic in the pion fields. Solitons in the first model are shown 
to deviate significantly from ansatze previously assumed in the liter¬ 
ature. The binding energies obtained in both models are lower than 
those obtained from the standard Skyrme model, and those obtained 
in the second model are close to the experimental values. 


1 Introduction 

The Skyrme model is a candidate model of nuclear physics in which nuclei 
appear as topological solitons. It can be derived from QCD as an effective 
theory valid in the limit of an inhnite number of colours [1]. It success¬ 
fully reproduces a number of qualitative features of nuclei, including quan¬ 
tum numbers of excited states [2] and the stability of the alpha-particle [3]. 
Quantitative features are generally predicted with reasonable accuracy, but 
with some notable exceptions. One class of quantities which are not accu¬ 
rately predicted within the standard Skyrme model is that of nuclear binding 
energies: typically, these are too large by an order of magnitude. 

In the last few years a number of variants of the Skyrme model have been 
proposed which aspire to rectify this problem. Among them are a holographic 
model that incorporates vector mesons [4], the near BPS sextic model [5], and 
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the lightly bound model [6]. This article contains a thorough investigation of 
binding energies in the latter two. These models have in common that they 
supplement Skyrme’s Lagrangian with terms involving the pion helds only. 

The Skyrme Lagrangian takes the form 

fc 

Here U{x) is an 5'f/(2)-valued held, = {d^U)U^ is its right-invariant 
current, and m-,, denote the pion decay constant and mass, and g is a. 
dimensionless parameter. This Lagrangian should be understood as dehning 
an effective held theory, and the ellipsis represents additional terms which 
are higher order either in derivatives or in the pion helds 1 — U. The baryon 
number of a held U{x) is dehned to be 

[ (^ijk^r{RiRjRk)d^x = ( 1 . 2 ) 

and is equal to the degree of the restriction of U to any spatial slice. 

The sextic model adds to Skyrme’s Lagrangian a term sextic in deriva¬ 
tives. The motivation for doing so is that a Skyrme model with only a sextic 
term and a potential term is BPS: energies are directly proportional to the 
baryon number, and hence binding energies are zero. Thus a Skyrme model 
that includes a potential and a sextic term with coefficients much larger than 
the other terms might reasonably be expected to have low binding energies. 

The lightly bound model does not include any terms of higher order in 
derivatives but instead includes an additional potential term proportional to 
Tr(l — f/)^. The model whose Lagrangian consists only of this term and the 
Skyrme term has an exact solution with baryon number 1 but has no solutions 
of higher baryon number which are stable to fission; nuclei are unbound 
in this model. One therefore expects the Skyrme model that includes this 
potential with a large coefficient to have low binding energies. Simulations 
of a two-dimensional toy model support this hypothesis [7]. 

In the sequel, numerical approximations to energy minima in both models 
will be presented for baryon numbers one to eight. These results are con¬ 
sistent with the hypotheses that low binding energies are attainable in both 
models. However, the sextic model has proved more challenging to simulate 
numerically and we have not been able to obtain reliable results for soliton 
energies in the parameter regime in which binding energies are expected to 
be comparable to those of real nuclei. The structures of energy minima ob¬ 
tained in the sextic model differ radically from those previously assumed in 
the literature [5, 8, 9, 10, 11]; thus a re-examination of the results obtained 


2 



in these papers seems to be warranted. The energy minima obtained in the 
lightly bound model differ from those in the standard Skyrme model in two 
respects: their symmetry groups are smaller; and they have the structure 
of ensembles of particles. The sextic model will be discussed in section 2 
and the lightly bound model in section 3; some conclusions will be drawn in 
section 4. 


2 The sextic model 

In this section we study a one-parameter family of models with energy 

Ee = j ^ - ^iRiRi + ^ (^^eijk"^i{RiRjRk) 

+ ( 2 . 1 ) 

e > 0 being the parameter. So, we have replaced the usual Skyrme term 
with a sextic term proportional to and modified the potential. The 
regime of interest is when e is positive and small. The choice of potential is 
motivated by two facts: skyrmions in the e = 0 model are compactons, which 
is numerically convenient (strongly localized solitons are less prone to hnite 
box effects), and this potential is massless. That is, the expansion of the 
potential about U = 1 has no quadratic terms, so the pions in this theory 
are massless. This is important because, were the pion mass non-zero, it 
would scale like and hence would be unphysically large in the regime of 
interest. 

It is convenient for our purposes to think of the Skyrme held as taking 
values in rather than SU (2). Thus we dehne ipo,... ,(p 3 so that 

U = ipol 2 + i'PlTi + i(p 2 T 2 + i(p3T3 (2.2) 

where ri, T 2 , T 3 are the Pauli matrices, and cp^ + ipf + ipl + ipl = 1. In terms 
of the held (p : ^ G the energy of interest is 

E, = ]- [ {e\d(p\‘^+ \(p*n\‘^+U{<pf)d^x, (2.3) 

^ Jr3 

where hi is the usual volume form on the unit 3-sphere and 

U:S^ U{ip) = l-ipo. (2.4) 
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This formulation of the model makes manifest the simple, but crucial, fact 
that Eq is invariant under volume-preserving diffeomorphisms of [5]: given 
any held ip and any diffeomorphism —)■ such that 

^jJ*{dxi A dx2 A dxa) = dxi A dx2 A dxs, (2.5) 

Eq{p o Ip) = Eo{p). It will be convenient to introduce the shorthand 

E^ = eE^^'^ + E^^^ + E^^\ ( 2 . 6 ) 

the superscript indicating the degree of the corresponding piece of the energy 
density as a polynomial in spatial partial derivatives. 

As is well known [12, 5, 13], in the case e = 0 this model is BPS: 


Eq > CuB, 


with equality if and only if 



p*uj = U{ip) dxi A dx2 A dxs. 


(2.7) 


( 2 . 8 ) 


or, equivalently, B = U. For our choice of potential, cu = Vol{S^) = 2n^ 
(since po = 0) and solutions of (2.8) have compact support, meaning 

0(a;) = V = (1,0,0,0), the vacuum value, for all x outside some bounded 
subset of Equation (2.8) can be interpreted as the condition that p is 
a volume-preserving map from the subset of on which (p ^ v to 
equipped with the deformed volume form Qlf = VL/U. From this we deduce 
that any B = 1 BPS skyrmion (solution of (2.8)) occupies a total volume of 


Voh = f % = dTT^ (2.9) 

and VoIb = BVoli, in obvious notation. There is a i? = 1 solution within 
the hedgehog ansatz. 


pnirn) = (cos/(r), sin/(r)n), (2-10) 

where the prohle function has /(O) = tt, /(oo) = 0, and satsihes the ODE 

df sin^ / 


2 =l-cos/. 
dr 

An implicit solution to this boundary value problem is given by 


( 2 . 11 ) 


-ir^ = /-7r + sin/. 


( 2 . 12 ) 
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There are many ways to construct charge B BPS skyrmions. One can, for 
example, superpose B charge 1 solutions (e.g. hedgehogs) with disjoint sup¬ 
port. Less trivially, one can precompose a charge 1 solution with a volume¬ 
preserving 5-fold covering map —)■ M^\M, for example 

ijB : ^ tIjB : (r, 6, 0) ^ (5-^/^, 0, 50), (2.13) 

in spherical polar coordinates. Previous phenomenological studies of BPS [5, 
8 , 9] and near-BPS [10, 11] skyrmions have all been based on BPS skyrmions 
of this type. Note that, by their very construction, these BPS skyrmions 
have spherically symmetric baryon density (they start with a hedgehog, then 
precompose with a map which preserves B). 

To simulate the model with e > 0 numerically, we put it on a cubic lattice 
of size with lattice spacing h, typical values being N ^ 15Iv^ and 
h = 0.05, so that the box length hN is, for each 5, roughly twice the radius 
of the BPS skyrmion given by (2.13). Dirichlet boundary conditions are 
imposed (so 99 = u on the boundary of the cube). We replace spatial partial 
derivatives by 4th order accurate difference operators to obtain a lattice 
approximant to E^. Starting at e = 1 with some initial choice of 99 , this lattice 
energy is then minimized using a gradient-based minimization algorithm (the 
quasi-Newton L-BFGS method, limited memory version). Having found a 
(perhaps only local) minimum of 5i, we then decrease e slightly and, using 
the El minimizer as initial data, minimize again. Iterating this process we 
construct a curve of local E^ minimizers, parametrized by e. To increase our 
chance of hnding the global minimizer of 5^, we repeat the process with a 
variety of initial guesses at e = 1 (generated by, inter alia, the rational map 
ansatz, and the product ansatz applied to sets of lower charge solutions), 
obtaining conjectured global minimizers for all charges from 1 to 8 . Using 
the same differencing scheme, we also construct a lattice approximant Bnum 
to 5. For all the results we will report we find that this is accurate (i.e. integer 
valued) to within 0.01 % but is always lower than 5. We take this as a rough 
measure of the expected accuracy of our energies. Finite box simulations 
generically underestimate energies because they exclude the contribution of 
the soliton tail. We compensate for this by reporting E = {B/Bnum)Ee. 

As well as monitoring the accuracy of 5, we also keep track of the Derrick 
scaling constraint, as in [14]. That is, any critical point of E^ should, by the 
Derrick scaling argument [15], satisfy the virial constraint 

D{e) := -^(35(®)(e) - eE^^\e) - 3E^^\e)) = 0. (2.14) 

Ee 

In a hnite box, large enough that boundary pressure is negligible, D{e) should 
be small and positive [14]. If D{e) becomes negative, or large and positive. 
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this indicates that our numerical solution is not reliable, and we discard 
it. In all the simulations reported here, we hud that D{e) remains small 
(less than 1.5%) for all e G [0.2,1] but then explodes shortly below e = 0.2 
when B > 1. Simultaneously to the bad violation of D{e) = 0, we hud 
that the numerical solutions develop spike-like singularities, indicating the 
development of spatial structures smaller than the resolution of our lattice. 
Once this occurs, we can no longer track the solution curve. We conclude, 
then, that the model is numerically inaccessible (at least with our methods) 
for 0 < e < 0.2. We will see shortly that, for B = 1 and B = 2, other 
methods can be used that allow us to get much closer to e = 0, but these 
rely on the enhanced symmetry of these cases. 

Let us denote by Emin{B,e) the energy of the global minimizer of of 
charge B (assuming that such exists). It is clear from the dehnition of E^ that 
this is a monotonically increasing function of e, and that E^iniB, 0) = 27r‘^B. 
It is bounded above by the total energy of any charge B BPS solution, for 
example, a superposition of B hedgehogs of disjoint support. Hence 

Emin{B, e) < B{271^ + eE<^^\^H)) ~ (19.74 + 72.36e)B, (2.15) 

so Emin{B, e) converges to 27r^i? as e —)■ 0. Figure 1 shows a plot of Emin{B, e) 
ioT B = 1,...,8, 0.2 < e < 1, as predicted by our numerical data. The 
energy value for H = 8 is taken from the minimum energy conhguration at 
each sampled value of e. Energy and baryon number data for e = 1, 0.2 can 
be found in appendix A, table 2. 

Figure 2 presents, for each degree from B = 1 to B = 8, a pair of energy 
minimizers, for e = 1 (left image) and e = 0.2 (right image). In each case, a 
surface of constant baryon density is shown, and the colouring of each point 
represents the orientation of ip = (pi, p 2 , Ps) there. Since the B = 1 mini¬ 
mizers are within the hedgehog ansatz, their pictures specify the colouring 
scheme: we radially project p onto the unit cube max{|(^i|, \p 2 \i Iv^sl} = 1 
and colour each face of the cube: where (±1, 0, 0) are yellow/cyan, (0, ±1, 0) 
are red/green and (0,0, ±1) are white/black. Except for B = 6, the e = 1 
minimizers are qualitatively similar to those of the usual Skyrme model with 
a pion mass: in both cases, the potential strongly penalizes helds where 0 
is close to —V over large regions, and hence, disfavours the hollow shell-like 
minimizers found in the Skyrme model without a potential, preferring the 
formation of chain-like structures. (The B = 6 skyrmion suggests that this 
effect becomes important at lower B for our choice of potential than for 
the standard pion mass potential.) Swapping the quartic Skyrme term for 
the sextic term does not seem to have made a qualitative difference to the 
skyrmions. A similar observation was made by Floratos and Piette [16], who 
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Figure 1: Minimum energies in the charge B sector for B = as a 

function of e. The dashed lines show the upper bound derived in equation 
(2.15). The dotted lines show the minimum energies for S = 1, 2 as a function 
of e when restricted to axial symmetry. 


found that the energy minimizers (for i? = 1 to 5) in the sextic model with 
no potential were qualitatively similar to those of the usual massless Skyrme 
model. We present two local minimizers for B = 8, labelled 8 a and 8b. At 
e = 1 it is clear that 8 a has lower energy (by 0.35%), but at e = 0.2 the 
energy difference is close to our numerical accuracy, and 8b is slightly lower. 

Comparing the e = 0.2 minimizers with their e = 1 counterparts, we find 
that as £ reduces, the skyrmions have a tendency to lose symmetry, the holes 
in the level sets of B tend to shrink, and B tends to become more uniform 
within the skyrmion core. The density of nuclear matter is roughly constant 
in real nuclei, so this effect is desirable, since it reproduces at a classical 
level a feature that may emerge only after quantisation in the conventional 
Skyrme model. It is along rays emerging through the holes, on which B = 0, 
that singularities begin to form as £ drops below 0.2. This is consistent with 
the hypothesis that, as £ —)■ 0 , the minimizers converge to low regularity 
solutions of the Bogomol’nyi equation (2.8): such a solution must develop 
singularities wherever B = 0 and 93 % u. A more detailed picture of this 
process can be seen in figure 3, which shows plots of B along a line through 
the centre of the B = 4 skyrmion, again for e = 1 and e = 0.2. This line 
pierces the cube diagonally through a pair of maximally separated corners. 
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Figure 2: Level sets of baryon density for skyrmions in the sextic model with 
e = 1 (left) and £ = 0.2 (right), degrees B = 1 to B = S. On each surface 
B = 0.5 max{i3(a;)}. 

Also plotted is the (square root) potential density U ior e = 0 . 2 , suggesting 
that (p is converging, as e —)■ 0 , to a continuous but not differentiable solution 
oiB = U. 

Except for B = 1 and B = 2, there is no evidence that ip is converging 
to a BPS skyrmion with axial symmetry. In the case B = 2, ip appears to 
converge to a BPS skyrmion with axial symmetry, but not one for which B is 
spherically symmetric. This casts doubt on the validity of phenomenological 
studies of BPS and near BPS skyrmions based on helds of the form ipn o V's 
(see equation (2.13)) [5, 8 , 9, 10, 11]. Such studies perform a rigid body 
quantization, and depend strongly on the isospin and spin inertia tensors of 
the BPS skyrmions used, 

= / d(po A d(pi A *(d(po A d(pj) + 0(e:), 

Jij{(p) = [ B{xf{r‘^5ij - XiXj)d^x + 0{e), (2.16) 

7e3 

respectively. For ipn ° 'ipB, J in particular is isotropic to leading order, 
J = const I 3 + 0{e), a very strong constraint which is not supported by 
the numerical data: see hgure 4. 

Since all the minimizers are axially symmetric for B = 1 and B = 2, the 



Figure 3: Baryon density on a diagonal line through the centre of the B = A 
skyrmion for e: = 1 and £ = 0.2. The dotted red curve shows W((^(x)) for 
£ = 0.2, where U is the square root of the potential (see equation 2.4). 



|.|^=l22 (numerics) 
I33 (numerics) 
Jii=J22 (numerics) 
J33 (numerics) 
I.|i=l22 (ansatz) 
I33 (ansatz) 
J.| .| =J22“*i33 (ansatz) 


Figure 4: Moments of inertia for the skyrmion with B = 2. The non-zero 
components of the leading terms in the e-expansion given in (2.16) are plotted 
for a range of values of e, together with the predictions from the ansatz 
O ^2- 

minimization problem can be reduced to two dimensions in these cases. 
We seek minimizers of the form 

(p{p,z,(j)) = {ni{p,z),n 2 {p,z)cosB(j),n 2 {p,z)smB(j),n 3 {p,z)) (2.17) 
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where {xi,X 2 ,x^) = (pcos0, psin0, 2:), and n is a map from the right half¬ 
plane to S'^ which maps the ; 2 -axis onto the great circle n 2 = 0. The energy 
of snch a map is 

E^{n) = 2 tt [ {e{\dpn\^ + \d^n\^ + ^nl) 

J{0,oo)xR I P 

'j 

+ ^{dpnid^n3 - d^nidpUsf -h (1 - ni)Hpdpd^. (2.18) 

This integral was minimised nsing a (zero temperatnre) simulated annealing 
algorithm. A logarithmic lattice was used for p with lattice spacings h = 0.05 
in Inp and ; 2 . The ranges of Inp and 2 ; were [-5,2] and [0,7.5], and partial 
derivatives were replaced by first order differences. In practice we found 
that minimizing a discretisation of the difference between the energy and its 
topological lower bound, rather than of the energy, improved the performance 
of the algorithm. 

As expected, the B = 1 minimizer is spherically symmetric, while the B = 
2 minimizer has deformed toroidal level sets of B. As e —)■ 0, the hole through 
the centre of these tori closes up, and the level sets develop singularities along 
the symmetry axis, see hgure 5. The reduction in dimension allows us to use 
a much hner grid, adapted to resolve the structures developing near the 
symmetry axis, so we can follow the solutions reliably down to much lower e 
than is possible with the fully three-dimensional, method. The virial quantity 
|D(e)| remained less than 1% for all minimisers that we found {D was not 
required to be positive in these simulations, as the use of a logarithmic grid 
for p effectively cuts out a neighbourhood of the xs-axis). The difference in 
energies of the axial and 3D simulation results were < 0.07% for B = 1 and 
< 0.18% for B = 2. 

Any static solution of the sextic model with e > 0 must satisfy an in¬ 
teresting property called restricted harmonicity [13]. If (p is a critical point 
of Ee, then E^ = + E^^^ + eE^^^ is critical with respect to all smooth 

variations of <p : —)■ and hence, in particular, with respect to all vari¬ 

ations of the form (p^ = tp o where tpt is a curve through the identity 
element in the group of volume preserving diffeomorphisms of Now, as 
already observed, and E^^'^ are invariant under volume preserving diffeo¬ 
morphisms, so it follows that itself must be critical with respect to such 
variations. This is precisely the condition for (p to be restricted harmonic, 
in the terminology introduced in [13], where it was shown that restricted 
harmonicity is equivalent to the condition that div(p*/i is an exact one-form. 
Here h denotes the Riemannian metric on the target space (in this case, S^). 
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Figure 5: Baryon density contours B = 1.8 in the plane Xi = 0 for the B = 2 
skyrmion at various values of 

Since = 0, this is equivalent to the third order nonlinear PDF 

d(div(p*h) = dkdi{diip ■ dj^p)dxk A dxj = 0. (2.19) 

For helds within the axially symmetric ansatz (2.17), this PDF reduces to 

p^{dpAn ■ d^n — d^An ■ dpu) + p^(9^n ■ d^n — dpU ■ dpd^n) 

= pOpU ■ dzTi — 2B‘^nidzni, ( 2 . 20 ) 

in which A = + 8“^. In hgure 6 we present evidence that our axially 

symmetric B = 2 numerical solutions are approximately resticted harmonic. 
The restricted harmonic map equation is a closed condition on the space of 
all helds (with respect to any sensible choice of topology thereon). Since 
minimizers of for e > 0 are restricted harmonic for all e > 0, it follows 
that, if the minimizers converge to some limit as £ —)■ 0, this limit should 
likewise be restricted harmonic. This is another reason to be sceptical of the 
axially symmetric BPS skyrmions used in [5, 8, 9, 10, 11]: it is known 

that, except for S = 1, none of these solutions are restricted harmonic [13], 
and that their failure to be so gets progressively worse as B increases. 
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Figure 6: Density plots for the left (a) and right (b) hand sides of the re¬ 
stricted harmonic equation (2.20) evaluated on the soliton with e = 1 and 
5 = 2, and their difference (c). The plots show that equation (2.20) is 
satished to within numerical error. 


Recall that the underlying motivation for this model is the hope that, for 
e small and positive, it may support skyrmions with classical binding energies 
close to real nuclear physics data. So is e = 0.2 small enough? Sadly, no. 
While 5 o.2 supports skyrmions with signihcantly lower binding energies than 
those found in the conventional (massless or massive) Skyrme model, they 
are still much larger than those found in nature: see hgure 11 in section 4. 
We can use the results of our more rehned axially symmetric simulations 
to predict roughly how small e would need to be to give classical binding 
energies in the right ballpark. Real nuclei have binding energies per nucleon 
of around 1%, so we should demand that 


BErnin (1, S') EjYiiyi (M 

.RF/mm(l) ^) 


0 . 01 . 


( 2 . 21 ) 


Imposing this for 5 = 2 yields e = 0.014, considerably beyond the reach 
of our main numerical scheme. (We use the rough hgure 1% rather than 
htting to the binding energy of the deuteron since the latter is anomalously 
small.) For e so small, it is reasonable to hope that (p will be close to a BPS 
skyrmion, that is, a minimizer of Eq = 5^®^ -|- The question is, which 

one? Except for 5 = 1,2 where enhanced symmetry allows a reduction in 
dimension, this question seems to be numerically intractable. Our results 
suggest that the minimizers converge as £ —)■ 0 to BPS skyrmions of low 
regularity, which typically have strings of conical singularities, and which 
preserve the property of restricted harmonicity. 
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3 The lightly bound model 

The lightly bound model is dehned by adding to Skyrme’s lagrangian (1.1) 
a term proportional to Tr(l — t/)^. After hxing units of energy and length, 
this model depends on two non-trivial parameters. In suitable units, the 
associated static energy functional takes the form 


E = 


(1 - a) ( --Tr(i?,i?,) + m2Tr(l - U) 


- -Tr([i?„i?,][i?„i?,]) + a(lTr(l - U)Y\d^x. (3.1) 


Here a G [0,1] and m := /F^^g) > 0 are dimensionless pa¬ 

rameters. In writing the energy in this way, we have implicitly chosen 
FT,/‘ig\/l — a as a unit of energy and — a/F t^q as a unit of length. The 

advantage of parametrising the model in this way is that soliton energies and 
sizes remain ronghly constant as a is varied, thus enabling simulations with 
a range of valnes of a to be performed on the same grid. In the limiting case 
q; = 1 the above model has a topological energy bonnd E > 87r^|i?| which 
can be satnrated only in the |i?| = 1 sector [6]. Solitons with |i?| > 1 are 
therefore nnstable to fission. The other extreme case a = 0 is the standard 
massive Skyrme model, in which binding energies are too large. 

Numerical approximations to minima of the lightly bonnd energy (3.1) 
with B between one and eight have been determined for a range of valnes of a 
and for m = 1. Fixing the valne of m ensnres that the Compton wavelength 
of the pion remains comparable to the size of a nncleus, as nnclear radii do 
not depend strongly on the value of a. 

To simulate the model for a given a nnmerically we used the same nu¬ 
merical scheme as in section 2 on a cubic lattice of size with = 301 
and lattice spacing h = 0.05. We obtained a lattice approximant to the en¬ 
ergy for a given valne of a, Fa for one or more (perhaps only local) energy 
minima in each charge sector. A wide range of initial conditions were used 
to assist in hnding the global energy minimimiser for each baryon nnmber 
considered. The accuracy checks, already described, were applied. Firstly, 
the baryon number was calculated using the formula (1.2) and evaluated us¬ 
ing the same scheme. The valne obtained was correct to within 0.001% of 
the integer valne. Secondly, the virial constraint D{a), similar to (2.14), was 
evalnated on energy minima (in which E^°'^ denotes the component of the 
energy involving a derivatives): 

D{a) := - E^^^ - = 0. (3.2) 

E 
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When evaluated on our numerically-determined minima this quantity was at 
worst 0.01%. 

For B between one and four the full range of a was explored. Two 
sets of simulations were performed: one in which the minima were initially 
determined for a = 0.95 and tracked as a decreased to 0 in steps of 0.05, and 
one in which the minima were initially determined for a = 0 and tracked as 
a increased to 0.95 in steps of 0.05. Both sets of simulations gave similar 
results for a close to 0 and 0.95. For small a the minima were qualitatively 
similar those of the Skyrme model and in particular shared their symmetry 
groups. For a near 0.95 energy densities of minima were localised at B points, 
which we interpret as the locations of the B nucleons within the nucleus. The 
transitions between the two types of soliton occured at different values of a 
according to whether a was increasing or decreasing. For a increasing the 
transitions occured when a G [0.4, 0.5], while for a decreasing they occured 
when a G [0.35, 0.4]. 


Configuration 

00 

V(^^)/=o 

BEPN (MeV) 

^y{r^)i=o (fm) 

1 

1.1082 

1.2457 

0.00 

0.8104 

2 

2.2119 

2.2261 

1.84 

1.4481 

3 

3.3129 

2.7627 

3.26 

1.7971 

4 

4.4111 

2.9310 

4.56 

1.9066 

5 

5.5132 

3.4593 

4.66 

2.2503 

6 a 

6.6128 

3.4426 

5.09 

2.4190 

6 b 

6.6130 

3.4426 

5.07 

2.2394 

7a 

7.7120 

3.5350 

5.45 

2.2995 

7b 

7.7127 

3.8888 

5.36 

2.5297 

7c 

7.7140 

4.1668 

5.21 

2.7105 

8 a 

8.8090 

4.0367 

5.95 

2.6259 

8 b 

8.8094 

3.9610 

5.91 

2.5766 

8 c 

8.8124 

4.1150 

5.59 

2.6768 

8 d 

8.8130 

4.1431 

5.52 

2.6951 

8 e 

8.8130 

4.1430 

5.52 

2.6950 


Table 1: Energies, isoscalar charge radii, and binding energies per nucleon of 
skyrmions in the lightly bound model with a = 0.95. 


Energy minima were examined more carefully for a G [0.9,1.0] and B 
between one and eight. This range of a is the one in which binding energies 
are comparable to experimental values. We found that a = 0.95 gave binding 
energies of appropriate magnitude relative to the proton mass; baryon density 
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Figure 7: Baryon density isosurfaces of skyrmions in the lightly bound model. 

isosurfaces for the minima that we found are displayed in hgure 7. Note 
that for B > 6 several local minima of the energy have been found. The 
isosurfaces have been coloured to indicate the orientation of the pion helds 
on the surface, as described in the previous section. The energies of these 
conhgurations are recorded in table 1. Also recorded in this table are their 
charge radii defined by the formulae 

(A.o = ^/^J|x-x„fB(x)d\, x„=^£xB(x)d»x. (3.3) 

The most striking feature of the baryon density isosurfaces is that they 
display a clear particle-like structure. This is in marked contrast with the 
traditional Skyrme model, in which individual nucleons are not discernible 
within a nucleus. The value of 0.1max{i3} used to create the isosurfaces in 
hgure 7 is signihcantly lower than those commonly used in the Skyrme model; 
if instead a larger value of B is used the particle-like structure is even more 
pronounced. The points at which B attains local maxima are close to the B 
points at which the held U takes the value —1. Unlike in the usual Skyrme 
model these preimage points are non-degenerate. The likely cause of this is 
the choice of potential in our model, which strongly disfavours U being close 
to —1. The standard potential Tr(l — U) has a similar but weaker ehect: 
with this potential —1 has two preimage points with four-fold degeneracy in 
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Figure 8: Subsets of the face-centred-cubic lattice corresponding to the iso¬ 
surfaces in hgure 7. 


the B = 8 soliton, whereas without it —1 has a single preimage point with 
eight-fold degeneracy. 

Another remarkable feature evident in hgure 7 is that, with one exception, 
all conhgurations resemble subsets of a face-centred cubic lattice. Within 
each of these lattice-like conhgurations, no more than four distinct orienta¬ 
tions for the constituent particles can be observed (corresponding to whether 
the particle is on a vertex of a cube, a horizontal face, or either of the two ori¬ 
entations for vertical faces). The precise subsets of the lattice are displayed 
in hgure 8, with the four orientations distinguished by four colours. The 
only soliton to which these comments do not apply is conhguration 8e. This 
conhguration resembles a pair of adjacent tetrahedra and cannot be realised 
as a subset of the face-centred cubic lattice; its eight constituent particles all 
have diherent orientations. 

In order to make a direct comparison between the lightly bound model 
and experimental nuclear physics we have calibrated our units of energy and 
length. The energy units were chosen so that the one-soliton had energy equal 
to the proton mass, 938.27MeV. The length units were chosen by htting the 
electric charge radii of the proton and neutron. If the proton and neutron 
are modelled as rigidly isospinning bodies centred at the origin then their 
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Figure 9: Binding energy per nucleon as a function of baryon number in the 
lightly bound model and from experimental nuclear physics. 


charge radii are given by 

(r2)p = f + 13 ) d^x, {r^)n = [ - X 3 ) d^x, (3.4) 

where X 3 denotes the isospin charge density, normalised such that Isd^x = 
We therefore identify 

^/{r‘^)i=o = {r^)p +{r‘^)n = VO.87832 - 0.1149 fm = 0.8103 fm. (3.5) 

These normalisation conventions force us to identify 

- ^ - = 10.72 MeV, ~ = 0.6505 fm. (3.6) 

Hence the following values are chosen for the parameters in the Skyrme la- 
grangian ( 1 . 1 ): 

= 29.1 MeV, = 303 MeV, g = 3.76. (3.7) 
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Figure 10: Isoscalar charge radii in the lightly bound model {a = 0.95), 
the standard Skyrme model (a = 0), and from experimental nuclear physics. 
The a = 0.95 curve passes through the points corresponding to global energy 
minima. 

With these units the pion decay constant is significantly below its experimen¬ 
tal value of 186MeV, and the pion mass is significantly above its experimental 
value of 138MeV. Most calibrations of the standard Skyrme model result in 
abnormally low values of the pion decay constant. Although this problem 
seems to be exacerbated in the lightly bound model, we expect for reasons 
to be explained below that the value of would increase in a quantised 
treatment of the model. A better value for the pion mass could be obtained 
by reducing the value of the parameter m in the model, but we have not 
explored this possibility. Reducing the value of m would probably reduce 
soliton energies, and thereby also increase the value of F.^. 

The binding energy per nucleon is plotted alongside experimental values 
in figure 9. The theoretical binding energy per nucleon is calculated from 
the formula Eb/B — i?i, in which Eb represents the global minimimum of 
the energy in the charge B sector. The corresponding experimental values 
are the binding energies for the lightest isotope at each mass number B. The 
theoretical curve closely follows the experimental curve for most values of B 
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in the range considered, except that the anomalies at -B = 4 and 8 are less 
pronounced in the theoretical curve than in the experimental curve. This 
graph represents a substantial improvement on the standard Skyrme model 
(see hgure 11 for a comparison). 

The skyrmion isoscalar charge radii are plotted alongside experimental 
values in hgure 10. For B even the experimental values shown are the charge 
radii of isotopes with mass number B and isospin 0. For B odd they are 
calculated as weighted averages of the states with isospin / = ±1/2 via the 
formula {r‘^)i=o = [{B + l)(r^) 7 =i /2 + {B — l){r‘^)j=_i/2]/2B. Like equation 
(3.5), this formula is motivated by identifying the I = ±1/2 states with 
solitons isorotating rigidly in opposite directions. 

Figure 10 shows that the theoretical and experimental charge radii are 
comparable. The most surprising feature of the experimental curve is the 
abnormally large radius at i? = 2 in comparison to B = 1. The lightly 
bound model does not fully capture this feature, but approximates it better 
than the standard Skyrme model, which is plotted alongside for comparison. 

4 Concluding remarks 

We have investigated the problem of obtaining low binding energies in two 
variants of the Skyrme model. Our results are summarised in hgure 11, 
alongside results obtained from a third “holographic” variant explored in 
[4], the standard Skyrme model, and the experimental binding energy curve. 
The sextic and holographic models are both capable, in principle, of produc¬ 
ing skyrmions with very low binding energies, but in practice they are yet 
to yield realistic binding energies in numerical simulations. For the sextic 
model, the problem is that the skyrmions start to develop singularities, and 
become numerically inaccessible, before the regime of very low binding energy 
is reached. For the holographic model, which is a Skyrme model coupled to 
an inhnite tower of vector mesons, only the lowest truncation level, including 
just the lightest vector meson and the lightest axial vector meson, has been 
treated numerically. One could, in principle, include more mesons from the 
tower, and this would presumably reduce the binding energies further. How¬ 
ever the resulting model would be formidably complicated, and hence difficult 
to solve numerically in practice. In contrast, the lightly bound model looks 
very promising: using standard numerical algorithms, binding energies have 
been obtained which are comparable with experimental values. This model 
clearly deserves further attention. 

The solitons obtained in the lightly bound model differ markedly from 
those obtained in the standard Skyrme model, in that they resemble a lat- 
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Figure 11: Binding energy per nucleon for B = 2,..., 8 in various variants 
of the Skyrme model; the massless Skyrme model with m = 0, the massive 
Skyrme model with m = 1, the sextic Skyrme model with e = 0.2, holo¬ 
graphic Skyrme model with vector mesons, the lightly bound Skyrme model 
with a = 0.95 and experimental data. Markers indicate values of available 
data. 

tice of nucleons. This result suggests that the problem of hnding solitons 
with higher baryon number could be simplihed by hrst investigating a point 
particle model, along the lines of [7]. 

Our analysis of the lightly bound model has been purely classical, but in 
the future this model should be quantised semi-classically. Including quan¬ 
tum corrections would no doubt modify some of our results. Most notably, 
the masses of quantised states would be higher than the classical values 
quoted here, with the one-soliton affected more severely than others. There¬ 
fore binding energies are likely to be reduced in a semi-classical treatment, 
and the value chosen for a would need to be reduced to compensate for this. 
Reducing the value of a would also lead to a more realistic value for the pion 
decay constant. 

Our attempts to obtain realistic binding energies within the sextic model 
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have been hampered by numerical difficulties associated with the formation 
of singularities in the held. Thus a different numerical scheme is called for 
if any progress is to be made with this model. Despite these difficulties, we 
have found clear evidence that solitons in this model do not resemble the 
ansatz proposed in [5, 8, 9, 10, 11]. In this light it would be worthwhile 
therefore to revisit some of the results obtained therein. 
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A Appendix: energies in the sextic model 


Conhguration 

e : 

= 1 


e = 

: 0.2 


El 

B 


Eo.2 

B 

1 

62.7543 

0.999991 


29.4172 

0.999995 

2 

116.696 

1.999981 


56.6878 

1.999972 

3 

169.146 

2.999970 


83.6968 

2.999764 

4 

220.424 

3.999964 


110.659 

3.999790 

5 

274.815 

4.999958 


138.065 

4.999874 

6 

327.445 

5.999950 


165.374 

5.999737 

7 

379.956 

6.999939 


193.238 

6.999792 

8 a 

433.719 

7.999931 


220.633 

7.999514 

8 b 

435.214 

7.999934 


219.884 

7.999600 


Table 2: Numerical values of energy and baryon number for skyrmions in 
hgure 2. is given by (2.3). 
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